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Abstract 

We consider a four-dimensional Af — 2 gauged supergravity coupled to matter fields. The 
model is obtained by a U(l) gauging of a charged hypermultiplet and therefore it is suitable 
for the study of holographic superconductivity. The potential has a topologically flat direction 
and the parameter running on this "moduli space" labels the new superconducting black holes. 
The phase diagram of the theory is studied and the zero temperature solutions are constructed. 
The model has rich dynamics. The retrograde condensate is just a special case in the new class 
of black holes. The calculation of the entanglement entropy makes manifest the properties of 
a generic solution and the superconductor at zero temperature is in a confined cohesive phase. 
The parameter running on the topologically flat direction is a marginal coupling in the dual 
field theory. We prove this statement by considering the way double trace deformations are 
treated in the AdS/CFT correspondence. Finally, we comment on a possible connection, in 
the context of gauge/gravity dualities, between the geometry of the scalar manifold in N = 2 
supergravity models and the space of marginal deformations of the dual field theory. 
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1 Introduction 

Over the past few years the AdS/CFT correspondence has been one of the most studied subjects in 
theoretical physics pQ. The duality has lead to the foundation of the holographic principle which relates 
a certain class of field theories at strong coupling to general relativity in asymptotically AdS spaces. The 
precise prescription for computing physical observables of the field theory through classical gravity makes 
the duality a powerful tool to deal with strong coupling dynamics [HE]- More generally, holographic 
techniques have been successfully applied on a vast number of problems, ranging from the quark-gluon- 
plasma to the most recent and challenging aspects of condensed matter physics [H[S] . 

The call for holography as guiding principle in condensed matter systems, comes from the lack of a 
field theory explanation for the properties of several many body systems. The physics underlying high-T c 
superconductors and strange metals seems to evade the well-established result that the Fermi liquid theory 
of fermions governs the low energy excitations of a Fermi surface, independently of the strength of the 
electron interaction at the lattice scale [6ll7] . On the other hand, this result remains an outcome of RG flow 
techniques in a regime in which perturbative quantum field theory can be applied. Holographic theories 
have been introduced as an alternative framework to overcome this issue. They naturally account for 
strong coupling physics and quite surprisingly, a recent work has shown that in a "realistic" holographic 
setting, some properties of the cuprates come for free. Thus, high-T c superconductors and strange metals 
may belong to universality classes that admit an holographic description [8UI2], 

The model presented in this paper belongs to the class of theories introduced as holographic supercon- 
ductors [T3HH]- The strongly coupled field theory is studied through its holographic dual and the field 
content of the Lagrangian, in the simplest case, consists of gravity, a complex scalar and a gauge field. 
Black hole geometries which support electric flux are dual to states of the field theory with finite density of 
matter. These geometries are generically unstable at low temperature and spontaneously develop a scalar 
hair which is responsible for a field theory operator to condense. Despite the simplicity of the model, the 
condensation seems to be a universal feature of charged gravitational backgrounds whenever a charged 
scalar field can be turned on. The only wrinkle in this story is given by the existence of the retrograde 
condensate, first pointed out in [TB]. This is an exotic and not thermodynamically favored phase of the 
theory, which has the odd characteristic of showing a condensate only for temperatures above a certain 
critical temperature. Nevertheless, it has been found in models coming from string theory. Motivated by 
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this last observation, we want to point out that the retrograde condensate has been discovered in models 
obtained by truncating a complicated supergravity theory to a simpler sector. Then, one might suspect 
that the simplicity of those truncations is hiding part of the story. In other words, we chase the idea that 
interesting physics comes together with the retrograde condensate. 

Our model is an J\f — 2 supergravity theory coupled to a charged hypermultiplet and it can be 
embedded in the four dimensional J\f — 8 supergravity. It includes two special solutions already discussed 
in the literature: a well-know holographic superconductor and a retrograde condensate. The main point 
of this paper is the construction of a class of hairy black holes whose condensates interpolate between the 
above mentioned special solutions. The new class of superconducting solutions present remarkable and 
novel dynamical aspects. In this sense, the idea that interesting physics was related (even if not directly) 
to the retrograde condensate is paid back. 

Af = 2 supergravity coupled to hypermultiplets is quite special. Hyperscalars are associated to 
coordinates on a quaternionic manifold and supersymmetry fixes completely the form of the Lagrangian. 
In this paper we analyze the dynamics of the two charged hyperscalars which parametrize the coset space 
SU(2, l)/{7(2). The scalar potential is non trivial: it has a local maximum at the origin, a saddle point 
and two tachyonic directions, one related to the saddle point and the other one related to the retrograde 
condensate; they both join at the boundary of the scalar manifold. The coset space is a two dimensional 
complex manifold with quaternionic structure but can also be described as an open ball with a non 
trivial embedding. Because of this topological property, a flat direction shows up in the scalar potential. 
An angle labels the set of AdS vacua associated with the flat direction and we call this parameter 9^. 
Then, the new class of hairy black holes is constructed by demanding that the superconducting solutions 
approach the asymptotically AdS vacuum at a fixed value of the parameter 6*00 . Tuning 6*00 along the flat 
direction changes the shape of the condensate. The well-know example of superconductor is found for 
9 = whereas increasing the parameter up to 9 = tt/2 turns the condensate into a retrograde condensate. 
In between, a cross-over of first order phase transitions appears. The phenomenology of these condensates 
is clarified by Figure l2~T1 and we refer to them as "interpolating solutions". The features of the condensates 
are not only qualitative like in a phenomenological theory: we quantify the properties of the interpolating 
solutions cooling the system down to zero temperature. 

In order to gain some intuition about the informations that the holographic principle provides when 
the theory is considered at zero temperature, we would like to stress the following important point. 
According to the Wilsonian renormalization group approach, low energy physics of a theory is obtained 
by integrating out the microscopic degrees of freedom down to the energy scale of interest. This idea 
is geometrically realized in theories with a gravity dual where the bulk radial coordinate represents the 
resolution scale of the dual field theory. In this case the short distances degrees of freedom of the field 
theory are associated to near boundary regions of the gravitational background whereas the low energy 
physics emerges from the shape of the geometry at small radial scales [17] . Holographic superconductors 
at zero temperature can be described in terms of RG flow dynamics. In the AdS/CFT correspondence we 
are implicitly assuming that the dual field theory lives on the conformal boundary of an asymptotically 
AdS space. In domain wall geometries between two AdS spaces, conformal invariance is also emergent at 
low energies |18j . In our extremal interpolating solutions a new scale appears in the gravity background 
and the geometry in the far IR collapses. We call this new scale C v . More precisely, below C n the bulk 
background has the shape of a cone all the way down to the tip. This characteristic brings to mind the 
gravity description of a confining theory as a cigar geometry [19] . Then, in the spirit of an effective field 
theory, we may wonder what kind of theory describes the low energy properties of the superconducting 
ground state dual to an interpolating solution. The holographic geometry suggests that this IR effective 
theory is a confining theory. To confirm this picture we study the behavior of the entanglement entropy. 

We recall that the entanglement entropy is a non local quantity which measures the amount of 
entanglement among the degrees of freedom of a ground state. It has a natural formulation in quantum 
many-body systems such as quantum spin chains [20)21 j . The system is characterized by a density matrix 
X and the entanglement entropy is given by the von Neumann entropy associated to the reduced density 
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matrix on a certain space domain £ , 

S = -Tr X £ logX £ , X £ = Tr £ c X . (1.1) 

The holographic dual of this calculation is understood as the analogous of the Bekenstein-Hawking en- 
tropy for a bulk surface which has d£ as its boundary [22] . The key point is that a transition similar to 
a confinement/deconfmement transition (see for example |23j ) takes place when the topology of the bulk 
surface changes from connected to disconnected. The result of our calculation shows that the new scale 
appearing in the zero temperature interpolating solution, sets the characteristic length £ of the confinc- 
ment/deconfinement transition. We use this transition as a measure of cohesion of the superconducting 
ground state. Indeed, in our zero temperature superconducting ground state, we expect that correlation 
lengths are finite and proportional to C v . This intuition agrees with the idea that the superconductor 
is better described as a cohesive state in a strong coupling regime [24] . Is this an indication that the 
superconductor will remain in a "glue" state at all temperatures? We believe that the answer is given 
by the order of the phase transition. A first order phase transition suggests that correlation lengths are 
finite also at the critical temperature. Our model is the first example of a field theory in which a tunable 
parameter, the angle #oo, accounts for all these phenomena. 

The confincmcnt/dcconfincment transition is the natural consequence of the geometry capping off 
below the new scale C v . We want to emphasize that the shape of the background in the far IR is the 
response (or the back-reaction) of the geometry to the presence of a UV mixing between the two charged 
scalars of the hypermultiplct. In the supergravity description this mixing is due to the non trivial 9oo. 
Interestingly, a similar mixing also relates the condensates which are dual to the two hyperscalars. This 
aspect is a remarkable outcome of the AdS/CFT scheme. Then, the field theory interpretation of the 
boundary data of our interpolating solutions gets interesting. Indeed, it is possible to prove that the 
parameter 6^ is the marginal coupling of a double trace deformation in the dual field theory. However, 
the deformation does not break conformal invariance and thus another mechanism must be responsible 
for the RG flow towards the IR confining phase. By considering the supergravity potential, it is evident 
that, in the classical picture, the mixing angle 6^ drives the hyperscalars along the tachyonic directions 
and therefore, towards the new IR geometry. We identify the field theory counterpart of this argument 
in the mixing between the two condensates: conformal invariance in the UV is broken and a relevant 
deformation initiates the RG flow. In particular, the same scale C v which measures the conical shape 
of the background, also measures how far an interpolating solution is from the IR conformal fixed point 
dual to the domain wall geometry. 

Several features play an important role in understanding the dynamics of the interpolating conden- 
sates. They are all tied up: the shape of the new extremal geometry, the confinement /deconfinement 
transition and the holographic RG flow interpretation. We believe that every aspect of our model has its 
roots in the properties of the scalar manifold and therefore, in the J\f = 2 theory under consideration. 
Ultimately, in its string theory origin. Unfolding these properties in the context of holographic supercon- 
ductivity increase the expectations that the AdS/CFT correspondence represents a fruitful approach to 
condensed matter physics. 

The organization of the paper is as follows. In section [2] we describe the main ingredients of the 
model and we define the holographic superconductor ansatz. In section [3] we discuss the details of the 
supergravity potential explaining the features of the topologically flat direction. We construct a more 
general class of hairy black holes presenting the numerical results for the finite temperature and the 
extremal solutions. We calculate the entanglement entropy to make manifest the physical properties 
of these new solutions. In section [4] we consider the field theory interpretation of our interpolating 
solutions. Using the holographic renormalization approach we explain how the parameter running on the 
topologically flat direction is mapped to a marginal coupling. Then, we give a detailed analysis of the RG 
flow of the theory. Finally, we relate the geometrical feature of the supergravity model to the existence 
of the marginal deformation in the dual field theory. 
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2 The model 



We consider the bosonic sector of AD N = 2 supergravity theory coupled to matter fields for a special value 
of the gauging parameter. The procedure to obtain the Lagrangian is illustrated in [2HJ12S3- The general 
structure of the theory includes: ny vector multiplets, rin hypermultiplets and the gravity multiplet. 
The scalars parametrize a manifold which is the direct product of a special Kahler manifold SK{ny) and 
a Quatcrnionic manifold Q(jih)- We analyze a simple model in this class of theories retaining only the 
gravity multiplet and one hypermultiplet. The quaternionic manifold is SU(2,\)/U(2) and the gauge 
symmetry is obtained by gauging a U(l) killing vector in the isotropy group SU{2) x U(l). The model 
is completely specified by the geometrical construction however, the gauging is not unique and the more 
general Lagrangian depends on the choice of a single gauging parameter. This parameter picks a U(l) 
direction in the isotropy group SU(2) xU(l). We consider the model obtained by gauging the 03 direction 
in the SU(2) subgroup. This particular model can be embedded in a consistent truncation of the J\f = 8 
supergravity [2 7) . 

The action of our model is given by, 

C = ( K- -F^Fi" - 2G mn (V»Q m {V~0n -V ) , (2.2) 

where F = dA is the field strength of the U(l) gauge field A. The quaternionic manifold SU(2, 1)/U(2) 
is a 2-dimensional complex manifold endowed with a metric G MN originating from a Kahler potential JC, 

K = -iln(l-|Ci| 2 -|C 2 | 2 ). ( 2 -3) 
In terms of the complex scalars £i and C2 the metric is explicitly written in the form, 

dG 2 = G MN dC M d( N , 

2 _ dCigr+jCggg (CirfCi + C 2 dC2)(CidCi + C2dC 2 ) r24l 

(i-iai 2 -ic 2 | 2 ) (i-icii 2 -ic 2 | 2 ) 2 ■ { ' } 

The covariant derivatives on the complex scalars are 

D^Ci = d^i + igA^t , D^2 = d^ 2 - igA^ , (2.5) 

and the potential is 



12 - 16(|Ci| 2 + | C2| 2 ) + 3(|Ci| 4 + |C 2 | 4 ) + io|Ci| 2 IC 2 | 2 

(1-IC1P-IC2I 2 ) 



As a consequence of the particular gauging of the rigid theory the scalars have opposite charges ±g. It's 
convenient to write the equations for £1 and £2 in the general form, 

(V M - igAJW - igA^Ci + (d M & - igA^X" - OV1C1 - , (2.7) 

(V M - igAnW - igA»)Q 2 + (d^2 - igA^X" - VV 2 Q 2 = , (2.8) 
where we have defined the quantities, 

2 



Xfj, — 



■(Ci^Ci+C 2 ^C 2 ) , (2.9) 



i-|Ci| 2 -|C 2 | 2 

Wl = v 4 r B l?i;-sy , (2.10) 



i-|Ci| 2 -|C 2 | 2 ' 

l-5|C 2 | 2 -3|Cip 

1-IC1P-IC2I 2 



w 2 = v 4 r 5 Sr^ 2 . (2.11) 



5 



The masses of the scalar fields can be read from the expressions of T>Vi and T>\>2- They are equal and 
are given by m 2 — —Ag 2 . 



2.1 Superconducting black holes: Ansatz 

The AdS/CFT correspondence in the present context is used as a tool to engineer states in the dual 
field theory with finite density of matter. We briefly recall how the dictionary between bulk fields 
and boundary operators works in our specific case. On the field theory side a finite density of matter is 
introduced by giving an expectation value (J*) to the charge density operator for a global U (1) symmetry. 
The current J M associated to the boundary global symmetry is mapped to a bulk Maxwell gauge field 

and the charge density operator is dual to the Maxwell gauge potential. The field theory is put at 
finite temperature by considering black holes geometries. The Hawking temperature of the black hole 
corresponds to the temperature of the field theory. 

The ansatz for the black hole is specified as follows. Requiring invariance under space-rotations and 
time-translations, the metric and the gauge field can be taken to be, 

t 2 dv 2 
ds 2 = -f(r)dt 2 + —{dx 2 + dy 2 + dz 2 ) + , A = $(r)dt . (2.12) 

L z J\r)h[ry 

Solutions of interest must be asymptotically AdS. As usual, the dual field theory lives in the UV 
region where r — > oo. In our model the AdS solution corresponds to the case in which all matter fields 
vanish. The corresponding metric components are f(r) = r 2 / L 2 and h(r) = 1, where we have chosen the 
normalization g 2 L 2 — 1/2. 

The holographic dual of a state with finite charge density is a black hole with electric flux at infinity. 
In this configuration $(r) has a non trivial profile. One candidate to describe this state is the Reissner- 
Nordstrom black hole, 

+ Kr) = l, S(r)=A,-£. (2.13) 

L z r 4r z r 

It represents an uncondensed phase because the charged scalar fields play no role in this background. The 
energy of the configuration is proportional to the mass M of the black hole whereas the charge density 
of the dual field theory is set by value of p. Looking for hairy black holes we assume the scalar fields to 
be functions only of the radial coordinate, 

Ci -> zi(r) C2 -> z 2 (r) e 1 ^ . (2.14) 

In this background the equations (|2.7[) and (|2.8I) are real and therefore the phases (pi and 4>2 are constants. 
We set them to zero for simplicity. Linearizing the equations of motion around the AdS vacuum we find 
the generic asymptotics, 

r 2 M 

f(r) = J2- — + - > ( 2 - 15 ) 
h(r) = 1 + ... , (2.16) 

*(r) = + , (2-17) 

r 

Zi(r) = + + . (2.18) 

The fall-off of the scalar fields is determined by the condition 

A(A - 3) = m 2 L 2 . (2.19) 
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The mass of the scalars is m 2 L 2 — —2 and therefore the solutions are A = 1, 2. In the general case, the 
boundary behavior for the scalar field can be written in the following form, 

£,(1) 

C=^A+7^, (2-20) 

where A is the larger of the two solutions of (I2.19j) . For masses in the range m 2 L 2 > m 2 BF + 1, the 
only normalizable mode is associated to the A. In the field theory language the scalar field is mapped 
to an operator of dimension A and the coefficient is the expectation value taken on the state dual 
to the gravitational background. The coefficient (of dimension 3 — A) is interpreted as a source for 
the dual operator and if the source is non vanishing then the field theory action is modified by a term 
proportional to 

SS oc / d 3 x 0^0 {2} . (2.21) 



This scheme goes under the name of "Standard Quantization". For masses in the range m 2 BF < m 2 L 2 < 
m 2 BF + 1, a slight modification of the action makes both modes normalizable and the role of source and 
condensate can be reversed [33]. In this case, the field theory operator has dimension 3 — A, the coefficient 
is interpreted as the condensate and is associated to the source. If the source is non vanishing 
then the field theory action is modified by a term similar to (12 .2 1[) . This reversed interpretation goes 
under the name of "Alternative Quantization". The mass m 2 L 2 = —2 belongs to the range in which the 
two types of quantizations are possible. 

The U(l) global symmetry of the field theory is spontaneoulsy broken by the condensate if the 
asymptotics of the scalar fields set the sources to zero 30 . Imposing regularity at the horizon allows a 
two-parameter family of solutions to the equations of motion. Once we require the spontaneous symmetry 
breaking condition, the superconducting black hole is specified only by the temperature. We leave further 
details to the Appendix. It is worthwhile to mention that we always work in the canonical ensemble where 
the charge density of the system is fixed. 



2.2 The old and the new 

The Lagrangian (12.2[) contains two superconducting black hole solutions already studied in the literature. 
We briefly review how they are embedded in our model (|2.2[) and what are their main features. 



• Model I, 

ft = tanh ~, C2 = 0, V = g 2 (sinh 4 (|) -4 (2 + coshry)) , (2.22) 
corresponds to [37] 

C = R - X -F 2 - ^{df]) 2 - 9 — sinh 2 n A 2 - V . (2.23) 

Choosing the alternative quantization scheme for Q\ a superconducting black hole is found for T c w 0.121, 
the phase transition is second order and the solution can be cooled down to zero temperature. It is the 
thermodynamically preferred phase. At T — the bulk geometry is nicely understood as a domain wall 
between two AdS^ regions with different radius. We will say more about this kind of geometry in the 
following sections. 



• Model II, 



Ci = -^tanh|, ( 2 = _L t anh|, V = -4g 2 (2 + cosh 77) , (2.24) 



corresponds to [2S] 

C = R-\F 2 - ^{dnf - 2g 2 Sinh 2 (|) A 2 V . (2.25) 
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The critical temperature of the hairy black hole solutions is the same of Model I, T c ss 0.121, but the 
superconducting condensate exists only for temperatures above T c . The solution is not thermodynamical 
preferred and the system always stays in the uncondensed phase. The shape of the curve suggests that 
the phase should be classify as a "retrograde condensate"!]]. 

Model I and Model II corresponds to two different scalar configurations however, they share a common 
feature, |£i| 2 + |^| 2 = tanh 2 (?7/2). This observation suggests that there is a larger class of black holes 
in our theorjU. In this new class of solutions Model I and Model II are only isolated points. In the 
next section we present a more sensible parametrization of the scalar manifold that nicely fits with our 
purpose. The situation is fascinating at least for one reason: it may clarify what kind of deformation 
relates Model I to the retrograde condensate. 



3 Interpolating solutions: Geometrical description 

The coset space SU(2, 1)/U(2) is topological^/ a ball and can be parametrized by a set of four real fields. 
Consider the following change of variables, 

Cl = T cos i e KvW>)/* f q 2 = t sin e -i( V -W2 . (3.26) 

It takes the metric (|2.4p into the form, 

d ° 2 = (T^F + W-^) {al + al) + W^?* 1 ' (3 ' 27) 

where g% are the standard left invariant one- forms |32) . 

&\ = cos ip d9 + sin ip sin 9 dip , 
&2 — — sin "0 d6 + cos ip sin 6 dtp , 
(73 = dip + cos 9 dip . 

The definition of the new variables is restricted to, r € [0, 1), 9 £ [0, 7r), (f € [0, 27r) and ip € [0, Air). We 
recall that the phases are always irrelevant; they do not appear in the potential and they are constants 
in the ansatz (|2.12[) . Defining r = tanh(fy/2) and A = A/y/2, we obtain the effective Lagrangian of the 
model 



£=^{K-\^-±smh*Q)(d9) 2 -l(dr ] f~lj(r 1 ,9)A,A^ V(r h 9)) . (3.28) 
The potential and the coupling J(rj, 9) in the field variables {rj, 9} are, 



P(t),9) = 5 2 (sinh 4 (|) cos 2 6> - 4 (2 + cosh 77)) , (3.29) 
J{rj,9) = 20 2 sinh 2 (-) (l + cos 2 6 sinh 2 (^jj . (3.30) 



t Related exotic aspects and a more accurate analysis of GL instabilities associated with a retrograde condensate can be 
found in |31| . 

'We thank Jorge Russo for this very important comment. 
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The equations of motion are 

*" + (M)*'- J <^>]!s= <" 3 » 

l" + (l + J + j)"' - jsMil «" + VM)^ - = (3.34) 

+ (I + 7 + £) +» th (I) + ^wiy (* J( "-'" m? - w a ' n "' e) ) - (M5) 



We discuss what values of (77, 9) are critical points of the potential V(r), 0). We need to solve the two 
algebraic equations, 

9 V V = -g 2 sinh ?7 (4 - cos 2 9 sinh 2 (|) ) = (3.36) 

<9 e 7? = -2 3 2 sin 26* sinh 4 (-) = (3.37) 

A simple solution is found for rj = 0. For this case d v V and dgV vanish automatically no matter what 
value of 9 £ [0, ir) we fix. Since the functions V(rj, 9) and J(tj, 9) are 7r/2-periodic it's sufficient to consider 
the following set, 

Me = {(17 = 0, 0) I e [0, tt/2]} . (3.38) 

The set A^e resembles a fiat direction but not in the usual sens^E Technically, the reason is that both 
scalars, Ci and £2, have non zero mass, m 2 L 2 = —2. What happens is that Aig appears as a flat direction 
only in the variables {77, 9}. Under our change of coordinates (I3.26[) . r\ — (or t = 0) implies £1 = £2 = 
and the whole Ale-set is mapped to the origin Ci = C2 = 0. Therefore, the interpretation of 9 is obtained 
by considering the ratio C2/C1 which gives a term proportional to tan(0/2). The choice of 9 fixes the 
direction of departure from the origin towards the boundary. This is a consequence of the topology of 
SU(2, 1)/U(2) that can be seen as the open ball in C 2 . Indeed, the coordinates (|3.26p . for fixed t^0, 
represent the Hopf fibration of the three-sphere. We refer to Aig as a topologically flat direction. 

It is useful to describe the behavior of the potential for a generic value of 0. In Figure [1] we show a 
density plot of V(j],0). The black dashed line is the set of points {(6,T)) | d v V = 0}. This line separates 
the plane (77, 0) in two regions according to the sign of d^V. In the colored region this derivative is 
negative and the potential decreases. In the white region the potential increases and it behaves like 
V ss exp 2?7 in the large r\ limit. In summary, Figured] shows that along the slices of constant 7^ 7r/2 the 
potential is bounded from below but, on the slice = n/2, the potential is negative definite and decreases 
like V ~ — exp 77. A second isolated critical point © exists for = and 77 = 2 arccosh-s/5 and according 
to the above analysis © is a saddle point. 

§In supersymmetric theories, as well as in supergravity theories, a classical flat direction in the potential is parametrized 
by a massless field whose expectation value sets the mass spectrum of the theory. 
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Figure 1: Density plot of the potential V ■ The potential gets steeper as the color becomes darker. The red dot is 
the saddle point 8 — 0, r\ — 2 arccosh\/5- For convenience the horizontal axis has been rescaled by a factor 0.5. 
The dashed black line is the set {(8,rj) | d v V = 0}. We anticipate the results of the section [3. 1\ The red lines, 
both dotted and solid, are the horizon values of the functions 6(r) and r/(r) when the superconducting black holes 
are constructed. In particular the value of Oca , which specifies the solution, can be read on the vertical axis: it 
coincides with 6(rh) when n(rh) ~ 0. On each line, the boundary conditions which define a superconducting black 
hole are satisfied. The figure shows that for the case 8(rh) 7^ the superconductor is driven towards Q[rh) = 7r/2. 



Independently of 77, the condition dgV = is satisfied by the two isolated values 8 = ir/2 and 9 = 0. 
These two particular values of 9 can be promoted to a constant bulk solution. Indeed, they correspond 
to Models I and II of the previous section: 

9 = Ci = tanh(r?/2), (2=0, (3.39) 



9 = n/2 Ci = 75 tanhfa/2), C2 - -\= tanhfa/2) , (3.40) 

The existence of the space Aig implies that 9 can be chosen arbitrarily when 77 vanishes. On the other 
hand 77 vanishes if and only if Ci =0 and £2 = 0. In a superconducting solution this last condition is 
satisfied at the boundary of the geometry, where the space is asymptotically AdS. Therefore it is natural 
to think that black hole solutions with fixed UV value of 9^ g A4g can be constructed. In particular, 
superconducting solutions with 0oo € Me interpolate between Model I and Model II. 

We explore this possibility by analyzing the asymptotic behavior of the variables {v(r),0(r)} in the 
background (|2.12j) . Expanding the equations of motion to lowest order in the 77 field, we find 



/' 1 6 $' 2 



= 0, (3.41) 



*"+(- + ^ ) *' = , (3.42) 



r 



h 



n" + (I + 4 + Pi rf - ]v 0- + ^77 + 4^ = 0, (3.43) 



r f h J 1 4 ' 4f 2 h 2 1 fh 

f h 



2 f h'\ 77' 

1 J ' " » °' + 2^-d' = . (3.44) 



r 



The equation for h(r) is trivial, h'(r) = 0. We suppose that the term 9' 2 is sub-leading. Thus the first 
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three equations decouple and are solved by the asymptotic expansion given in (|2.15[) - (|2.18[) . In particular, 



r)(r) 



0« , of 



+ 



(3.45) 



Substituting into (|3.44p we find that the choice of quantization scheme for rj(r) crucially affects the 
asymptotics of 0(r). The equation, 



has the desired solution only if Or, = 0. Only in this case we find the admissible asymptotic behavior, 



As 0' 2 is of order 1/r 4 we find that the assumption we made was correct. At this stage, we can think of £ 
as the condensate relative to 9^. In the next section we construct superconducting black holes with the 
following boundary conditions: 9^ £ Jv[ s and = 0. The coefficient 0|j = O v define the condensate. 
These are the interpolating solutions that we want to construct. Our solutions offer a rich (and quite 
unexpected) framework where all the physical properties associated to them have a simple explanation. 

3.1 A = 1 Condensation: Numerics at finite temperature 

We study the full system of equations (|3.3ip - (|3.35l) with the asymptotic conditions defined in the previous 
section. The condensate O v — O^T^Ooo) is associated to an operator of dimension A = 1 and it is 
function of the temperature T and of the angle Ooq. Figure |2~T1 shows the temperature dependence of the 
condensate when the value of 6^ is varied in the range [0, tt/2]. 

All the curves originate from the same branch at T ps 0.121. This fact is easily understood consid- 
ering the linearized approximation around the Reissner- Nordstrom black hole [33]. Indeed the critical 
temperature is only determined by the quadratic terms in the equation of the charged scalar r)(r). In 
the present case all the terms depending on the field 9(r) are suppressed and therefore only the mass, 
m 2 L 2 = —2, and charge, qL = 1/2, enter in the analysis. It follows that the temperature at which a 
small superconducting black holes exists does not depend on 9^,. 

Important changes happen when the value of 9^, is tuned up to the value 7r/2. As we expected, the 
shape of the curve O v (T, Oca) gets closer to the retrograde condensate, but for intermediate values, for 
example 9^ = 1, something interesting happens. In this case the shape of the curve brings to mind 
the case of a first order phase transition [51] , Calculating the free energy of the superconducting black 
holes and comparing the result with the Reissner-Nordstrom case, we find that the order of the phase 
transition depends on 9^ . The situation is clarified by Figure 12.21 In the white region the uncondensed 
solution is the only allowed configuration. When the temperature is cooled, the superconducting phase 
is thermodynamically favored. In the blue region the phase transition takes place at T c sec ps 0.121 and 
it is second order. It means that the system enters the superconducting phase in a continuos way. This 
mean field behavior is a consequence of the large N limit. In the yellow region and for temperatures 
above T c sec the phase transition is first order. The critical temperature of the first order phase transition 
can be read from the contour of the yellow region: it increases along the dashed red line. In this case the 
curve of the condensate has two branches, the lower branch is not thermodynamically favored and when 
the temperature approaches the critical temperature the system jumps into the superconducting phase. 
The free energy shows a discontinuity similar to the cases studied in [3"4"] . 

When the phase transition is still second order something new happens, even if the system has entered 
the superconducting phase, the free energy shows a discontinuity at some temperature smaller than T£ ec . 
The new behavior starts approximately at 9 crlt ps 0.95 and the discontinuity is of first order type. The 
discontinuity is physically interpreted as a change in the nature of the degrees of freedom that are 
responsible for the phase transition. Indeed, a pictorial feature of the curves 0^(T, Oqo) when 9^ is not 




(3.46) 



6{r) = 0oo + - + ■ ■ • • 



(3.47) 
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Figure 2.1 



Figure 2.2 




Temperature 



Figure 2: Numerical results for finite temperature superconducting black holes. The normalization has been fixed 
to p — 1. On the hand left side the condensate O v (T, Boo) is studied by letting the value of 6 X vary. The black 
dots are the values of the condensates extrapolated from our new zero temperature solutions. From bottom to top 
the curves correspond to Oao = 0,0.1,0.25,0.5,0.75,1. In the cases 9oo ~ and Oca = 0.1 the zero temperature 
values of the condensates differ by a factor of about 0.01. The rightmost red curve is the retrograde condensation 
which corresponds to 9oo = n/2. Strictly speaking the curves should be represented in a three dimensional space 
according to their value of 9 X . For simplicity we have collected the results in a single picture. On the right hand 
side the figure shows the phase diagram of the system. 



small is the presence of a plateau that starts at some intermediate value of the temperature and describes 
the remaining part of the curve as the zero temperature is reached. The extension of the plateau grows 
with and when the phase transition is first order it basically represents the entire curve. In the next 
sections we give more details about these phenomena by calculating the entanglement entropy of the 
superconducting ground state. Before, we discuss the existence of the zero temperature solutions. 



3.2 Zero temperature solutions 

Figure 12.11 indicates that for several values of 9oo £ M.$ the numerical black holes converge to a zero 
temperature solution. The case 9^ = is well studied and for completeness we review the argument. 

In a superconducting solution the value of the rj(rh) generically runs along the potential. Small values 
of r]{rh) correspond to temperatures in a neighborhood to the left of the critical temperature. In the case 
0oo = 0, the field 9(r) vanishes in the bulk and thus the value of rj(rh) is tied to the slice {77, V(rj, 0)}. 
When the temperature is cooled, rj(rh) hits the critical point &. The situation is graphically realized by 
the solid red line in Figure [TJ At the critical point there is an emergent A0IS4 space with radius L 2 = 3/7 
and the solution of the equations of motion is, 



7 2 

3 f 



h{r) = 1, i](r) — 2 arccosh-\/5, <I>(r) = 



(3.48) 



Because the gauge field must carry flux at the boundary, we have to move from (|3.48l) exciting irrel- 
evant perturbations, i.e sub- leading modes of the fields $(r) and rj(r). On this background the small 
perturbations, Srj(r) and <5$(r), are governed by the equations 



5$" + _«$$' 



and the solutions we need are, 



60 5^> 



. „ 4. , 60 Sr] 
r 7 r z 



, 



r)(r) 



T)IR r 



$(r) = §ir r + 



with 



(3.49) 
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The extremal solution is obtained integrating the full system of equations up to the UV boundary where 

(2) 

the shooting method is applied to match the condition O r) = 0. The background resembles a domain 



wall between the emergent AdS^ and the boundary AdSi |35U36| . Sometimes we refer to this solution as 
the conformal domain wall. 

For the cases Qoo 7^ one could be suspicious about the zero temperature limit. We have already 
observed that the potential has a runaway direction from the saddle point 6 towards the value 9 = 
tt/2. Exactly for 9{r) — tt/2 the condensation is retrograde and one might think that all the curves 
corresponding to 9^ ^ have to be retrograde at some point. Indeed Figure [T] shows that for fixed 
< 6>oo < tt/2 the numerics of the solutions when T — > converges to 8(rh) — > tt/2. In the same 
regime Tj(rh) grows and $'(r^) takes small values. On the other hand, if the extremal limit of the finite 
temperature black holes exists, we should look for a solution of the equations of motion with the following 
behavior: 9(r) — tt/2, $(r) = and divergent scalar field rj(r). Surprisingly, an analytic solution can be 
found and it is given by, 

dr 2 C it 

ds 2 = r 2 (-dt 2 + dx 2 ) + — — , ij(r) = 2 arcsinh -3. , 9(r) = - , $(r) = . (3.50) 

v t 2 



From ds 2 we read the definitions of the metric functions, f(r) = r 2 and h(r) = + C 2 /r 2 . We are now 

in the position to construct the extremal black hole. We look for irrelevant perturbations 68 (r) and 5$(r) 
on the background (I3.50p . Linearizing the equations for $(r) and 9(r) we find that the perturbations are 
governed by the equations, 



6$" + ~6& - ^<5$ = , 89" + -69' - -^58 = , (3.51) 



and the solutions are, 



8$(r) =CPV+^- , 66{r) = C^ 1) r+^— . (3.52) 



Imposing = and Cg — 0, we integrate the full system of equations, from small radial scales up 

to the UV boundary. There are three parameters left, {C v , Cg, C$}. We use the radial scaling to set C v 
to unity and the shooting method technique to match the boundary conditions. In particular we want to 
fix 9oo in the range (0, tt/2) and we need to impose — 0. From these new zero temperature solutions 



we extrapolate the values of the condensates. Figure 12.11 shows the result, the zero temperature value 
of the condensate, for fixed #oo, agrees with the limit O n (T — > 0,#oo) taken from the finite temperature 
black holes. 

In the approximation C n — > the metric in (|3.50l) matches with AdS^ and the divergent behavior of 
the scalar field is turned off. One might wonder if in that limit the expansion (I3.48[) is somehow recovered. 
We want to point out that this is not the case: for arbitrary small values of 0^ > the limit 9(r —> 0) 
is always tt/2. The correct interpretation of the statement is to think about the IR values, 8 = and 
8 = tt/2, as "order parameters" that distinguish two different phases of the theorjQ. This consideration 
implies that the IR physics described by the conformal domain wall and the new zero temperature solution 
is different. Indeed, a new scale appears in the solution: C v . The first geometrical feature related to 
C n is evident: in the region r -C C v the metric looks like a cone in which the transverse space M is the 
Minkowsky space, 

ds 2 w dr 2 + C 2 r 2 {-dt 2 + dx 2 ) . (3.53) 

The cone structure can be seen explicitly through a formal construction, 

• For < C v < 1 we define cos a = C n and we consider the following embedding. The curve 
z/r = ±tana with r > 0, defines a "surface of revolution" in the 5-dimensional space given by 
ds 2 = dz 2 + dr 2 + r 2 dM 2 . The resulting 4-dimensional metric is ds 2 = dr 2 



"We thank R. Myers for a discussion which led to this line of thought. 
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Figure 3.1 



Figure 3.2 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.0 0.5 1.0 1.5 



6(oo) e(<x>) 

Figure 3: From right to left, the functions C^((?oo) and O v {T — 0,800). The zero temperature solutions are 
normalized taking p = 1. According to our numerical precision, the best fit of the curve O v (T = 0, #00), for small 
values of Boo, is [O v — 1.265) « 0.55 #S' 75 ' ) . It has been obtained by considering a log-log plot of the curve. The 
values of O v (T — 0,#oo) not available from the numerics of the zero temperature solutions can be obtained taking 
the limit of the corresponding finite temperature black holes. Even in this case the fit agrees with the numerical 
data. On the right hand side, the dashed black line highlights the value C v = 1. As mentioned in this section \37B 
the value C v = 1 distinguishes between the two 5d embeddings of the IR cone geometry. 



• For C n > 1 we define cosh a = C n and we consider a similar embedding. The curve z/r = ± tanha 
with r > 0, defines a "surface of revolution" in the 5-dimensional space given by ds 2 = —dz 2 + 
dr 2 + r 2 dM 2 . The resulting 4-dimensional metric is ds 2 — dr 2 + r 2 cosh 2 a dM 2 . 

The construction is formal in the sense that, unless one of the transverse direction in M is compact, 
there are no angles that can be used to rotate the curve z = z(r). At r = 0, the bulk caps off with 
a "good" conical singularity. Indeed, according to the usual criteria, the on-shell potential is bounded 
from above [37]. We can check the statement looking at the value of the potential at origin r = 0. In 
this case 9(0) = tt/2 and V{r)(r),ir/2) = -2(2 + cosh 77) = -2(3 + 2C 2 /r 2 ) goes to -00 in the limit 
r — > 0. On the other hand, for intermediate scales C, < r < r, where f is such that deviations of the 
bulk solution from (|3. 501) are important, the metric still looks like AdS4,. A quantity that may be able to 
distinguish the differences between the conformal domain wall and the new zero temperature solution is 
the entanglement entropy We dedicate the next section to this calculation. 

Figure |31 shows the relations between C v , O v (T — 0) and O^. We already know two limits of these 
curves. The first one is obtained for 9^ —> 0. In this case the solution is close to the conformal domain 
wall, we can check that C v is converging to zero whereas the condensate is approaching the expected 
value O v (T = 0,9^ = 0) w 1.265. The second limit is 9^ —> ir/2. We know that the condensation is 
retrograde and according to the numerical data the curve has no turning point towards zero temperature. 
As a matter of fact, both values of C v and O v (T — 0) become singular meaning that the extremal solution 
blows up. This behavior agrees with the expectation that the curve of the condensate exists only for 
temperatures above T c . 

The last observation about the properties of the solution regards the electric flux. The electric flux 
T(r) is given by the expression, F{r) = r 2 h(r)<&' (r) . In the region r <C C v it behaves like F(r) s» C^C$r 
and therefore, in the limit r — > 0, the flux vanishes. From this observation we learn that the density 
charge p of the field theory originates only from the charged matter in the extra coordinate. Indeed, 
integrating once the equation p. 331) we obtain the Gauss law for the superconducting solution. It re- 
lates the U(l) charge p, the electric flux at the origin and the integral of the coupling J(r],9) across the 
bulk. Since the flux vanishes at r — 0, the field theory charge p originates only from the integral of J(ry, 9). 
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Figure 4.1 Figure 4.2 




Figure 4: Cartoon of the bulk surface 7£ used in the calculation of the entanglement entropy. On the left hand 
side the surface is connected and the configuration has a turning point located at z = z+ . On right hand side the 
surface is disconnected, the two planes are located at fixed x-coordinate and they are extended in the z direction, 
from the boundary up to z = z*. 



Observation. The zero temperature solutions that we have constructed answer a puzzle, pointed out 
in [35] , about the retrograde condensate. We briefly review the argument and we refer to [38, for further 
details. The authors study families of charged solitons, defined by various boundary conditions, in the 
context of J\f — 8 supergravity. These solutions can usually be parametrized by their mass. In the cases 
in which the mass is unbounded and the soliton can get arbitrarily large, a planar limit exists. This 
limit corresponds to the zero temperature solution of the holographic superconductor which belongs to 
the specific theory under consideration and is defined by the same boundary conditions of the solitons. 
Model II was analyzed and it was found that the planar limit of A — 1 charged solitons is given by (|3.50[) . 
but the zero temperature superconductor doesn't exist because of the retrograde condensate. Now we 
understand that Model II is just a special case of 9^ £ M.q. Nevertheless, as we have seen in this section, 
the solution p.50p represents the IR limit of a large class of extremal black holes with each of them being 
labeled by 8^. 



3.3 The Entanglement entropy and the confined cohesive phase 

The entanglement entropy of a region £ in the boundary field theory is computed holographically accord- 
ing to the proposal of [35] ■ It is calculated by considering the area of bulk surfaces j£ whose boundary 
is given by d£. According to the proposal, the entanglement entropy is the minimal area, 

S £ = 27rAre 2 a(7£) • (3.54) 

The strategy is to rewrite the problem as a variational problem. The solution of the equations of motion 
provides the profile of the minimal area surface and allows one to compute the entanglement entropy of 
the configuration. It's convenient to use the variable z = 1/r writing the spatial part of the metric in the 
form 

L 2 1 

Spatial = ~£ {dx 2 + dy 2 + U (z)dz 2 ) , U{z)= z2f ^ h2 ^ • (3-55) 

The UV boundary is located at z = 0. We are interested in boundary surfaces with strip shape, 



£ = {(x,v)\- l ±<x< l -Z,0<v<l v }. (3.56) 



We parametrize the bulk surface j£ choosing coordinate y and x — x(z) as schematically depicted in 
Figure 14.11 The area functional is given by, 
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Area( 7f ) = 2L 2 l y J — ^U{z) + x'{zf . (3.57) 

In the formula is the turning point of the configuration whereas e is an UV regulator. The variational 
problem for x(z) has a conservation law that allows to eliminate x'{z) from the integral f|3 . 57|) . The 
conservation law is, 

_ = 1 f 3 5 8 ) 

Z 2 y/U{z)+X , (z)* Zl' 

We solve for x' (z) and we use the result to write the area functional in the form, 




Area( 7£ ) = 2L 2 l v / - J \ . (3.59) 



Integrating the conservation law we find the relation between l x and z+, 

2" 




(3.60) 



The entanglement entropy is divergent in the limit e — > 0. The origin of the divergence is easily under- 
stood: it corresponds to the integration of the short distances degrees of freedom and geometrically is 
due to the fact that the minimal surface reaches all the way to the boundary. This leading divergent 
term is the "boundary area law" of the entanglement entropy. Subtracting the divergence, we can write 
Sg in terms of the finite quantity s defined by, 

In the following the curve s — s(l x ) is computed using numerical calculations. Formula p.61[) is propor- 
tional to the ratio L 2 /n 2 , where L is the AdS radius and 1/k 2 is the gravitational constant of the 4D 
model. According to the gauge/gravity correspondence L 2 / k 2 is a function of the numbers of colors of 
the dual gauge theory. The present model is not based on a string construction thus, the precise relation 
is not fixed. However, in the large N limit the dependence L 2 /k 2 ~ A 3 / 2 is to be expected [40] . The 
overall coefficient is not determined but the N 3 / 2 dependence is a robust feature. 



Finite Temperature. We begin by analyzing the results at finite temperature. A general discussion 
about this case has been presented in [41] and here we recognize basically the same features. A distinction 
between first order and second order phase transition in the condensate phase diagram is necessary. 

We first consider the case in which the phase transition is second order. In Figure [5] we show the 
plots of s = s(l x ) for the two values, 6*00 = 0.1 and 6^ — 0.5. When the temperature is close to the 
critical temperature, the large lengths behavior of the entropy shows a linear dependence. Lowering the 
temperature, the slope drops to zero and the entropy approaches a constant value. If the region E has a 
relatively small size then is close to the boundary and thus the pure AdS^ result l/l x is recovered. 
In section 13.11 we observed that at fixed low temperature, depending on the value of #oo , the condensate 
Or)(T, 9oo) enters the plateau region. In this cases the linear behavior s(l x ) ~ l x drops very fast if 6^ is 
not small. An example is shown in Figure IS~21 

When the phase transition becomes first order, the entanglement entropy is multi-valued at some 
length Ik a kink appears. The swallowtail curve, typical of this cases, is displayed in Figure 16.11 For 
I > Ik the slope of the curve s(l x ) suddenly changes and the entropy saturates to a fixed value. 
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Figure 5: Numerical plots of the entanglement entropy as function of the length l x . On the left hand side the 
case 800 = 0.5 is displayed for different temperatures. From bottom to top T ~ 0.04,0.095,0.11. On the right 
hand side, we show a comparison between the behavior of the entropy as function of 9^ for a fixed temperature 
T ss 0.04. The dashed line corresponds to 9^ = 0.1 whereas the solid blue line is 9^ — 0.5. The slope for the 
curve corresponding to 9^ = 0.5 is approximately constant for large l x . 



Zero Temperature. We now turn to the zero temperature solutions. For the conformal domain 
wall the picture does not show any substantial novelty with respect to the general analytical arguments 
reviewed in |42U43j . We consider directly our new numerical results. We repeat a simple scaling argument 
that captures the main feature of the zero temperature solution when 9^ ^ |44j . 
If z 1 the function U(z) is given by the radial component of the metric in (I3.50[) , 

U(z) = (1 + C^z 2 )- 1 ps (C^z 2 )" 1 z>l/C r , region Ri, (3.62) 
ps 1 z < z < l/C v region R 2 , (3.63) 

where z — 1/r. If the length l x is probing the intermediate region R2 then we expect the known behavior, 
however deep in the IR something new happens. We can estimate the integral (|3.60[) by considering 
(I3.62[) and the change of variables z — > z/z*, 

I f 1 z 2 ( 1 \ 1 

-W = z * I dz U(zz*) 1/2 ~ const z* I — — I ~ const — . (3.64) 

* JO vl-Z 4 yZ+Urj/ Crj 

A maximum length l max exists and the above relation shows that it is proportional to C^ 1 . According to 
the interpretation of the radial coordinate as energy scale, Z max and C v are correctly related. It might seem 
surprising that in the limit of large z+ the length l x remains tied to the value V^ ax . However, it does not 
mean that there are no configurations with l x > Z™ a:c . It is important to keep in mind that the previous 
calculation considered only smooth and connected surfaces 75, but in addition we have disconnected 
configurations. This class of surfaces consist of two disconnected planes located at x = ±l x /2 that are 
extended in the z direction for a length equal to z*. The entanglement entropy as a function of the length 
l x is constant and only depends on z+ through the following formula, 

4ttL 2 , f z * dz r—— AttL 2 , / 1\ . . 

S z = — l y] Q ^VW) = ^rly(s+-) ■ (3-65) 

Figure IB~2l shows that when the size of £ is stretched up to Z™ ax , there is a transition from the con- 
nected to the disconnected surface. In the limit z* 3> the contribution to the entanglement 
entropy coming from the disconnected configuration converges to the value reached by the connected 
surface when l x = l max . This transition is similar to a confinement /deconfincment transition [45 . The 
rate of change of the entanglement entropy with the length l x , diS, makes clear the connection with 
the confinement/deconfinement transition. The features of the transition are opposite from that of an 
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Figure 6.1 



Figure 6.2 




Figure 6: Behavior of the entanglement entropy for several cases. On the left hand side we show a zoom of the 
swallowtail curve which characterizes a first order phase transition, 9 X = 1 andT w 0.082. On the right hand side 
we present the result for the new zero temperature solutions. Along the solid lines the entropy has been calculated 
making use of the connected surface 'ye ■ Dashed lines correspond to the contribution of the infinitely long discon- 
nected surfaces. The maximum length l max is indicated by the black dots. From top to bottom the various cases 
correspond 9^ — 0.5, 1, 1.4. At l max the solid lines merge with the dashed lines and the confinement/ deconfinement 
transition takes place. Interestingly, for 9^ = 1, Ik is of the same order as l max . 



Hawking-Page transition [23]. For connected configurations diS ~ iV 3 / 2 and the degrees of freedoms that 
form the ground state are in a deconfined phase. When the topology of the bulk surface changes to the 
disconnected configuration, the degrees of freedom living in £ are not correlated with the ones of £ c . 
This situation characterizes the large length scales and indeed the entanglement entropy is a constant for 
l x > l" lax : for disconnected configuration diS vanishes. In this sense we think of l max as a sort of cohesion 
length which the entanglement entropy is able to probe. According to the classification introduced in [24], 
a confined cohesive phase is emerging in the IR of our model. 

At this point we can give a more exhaustive answer to the question that we posed in the introduc- 
tion: what is the effective field theory that describes the low energy properties of our zero temperature 
superconductors? In the case Soo ^ 0, this field theory has the characteristics of a confining theory. The 
geometry of the interpolating backgrounds is effectively cut off at the scale C v and a log term appears 
in the IR expansion of rj(r), i.e. rj(r) s» log(C 2 /V 2 ) for r <C C v . The confined cohesive phase which 
emerges in the IR of our solutions is the most fascinating and intriguing novelty: the back-reaction to the 
presence of a mixing between the two charged hyperscalars is responsible for the confining IR geometry. 
On the other hand, when 8^ = 0, the end point of the domain wall is dual to a conformal field theory. 
Concerning this, it is interesting to observe that l x ~ 1/C V goes to zero in the limit 8^ —> 0. 

The calculation of this section helps to understand what characteristic length can be associated to the 
superconducting state in order to distinguish between first or second order phase transition. We consider 
the case 8^ = 1 which is particularly instructive. The phase diagram of O v (T, 1) shows a clear first 
order phase transition. It means that correlation lengths are finite when the temperature approaches 
the critical temperature. Because there are only two scales in our gravitational solutions, CV, and T c , all 
other scales will be functions of C n and T c . Then, we can associate the (zero temperature) value of l max 
to this first order phase transition. 

For a generic 8oo we need a more careful analysis. Indeed, we know that the maximum length l max 
exists for each value of #oo but the phase transition is not always first order. In particular, for small 
values of #oo the scale C n <C 1 and therefore Z max ^> 1. In these cases the phase transition is second 
order. We can make a qualitative comment about this phenomenon recalling the discussion of section 
13.11 When 8^ increases, the curve of the condensate develops a plateau (see Figure l2~Tj) . This plateau 
grows until the phase transition becomes first order. At the same time l™ ax decreases towards a critical 
value related to 8 cnt . This critical value of the length l™ ax is the one that we associate to a "strong" 
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cohesion of the ground state. 

In summary, by calculating the entanglement entropy we have unfolded the interesting features codihcd 
in our extremal metric p.50p . The scale C v sets the characteristic length of the confined cohesive phase. 
The appearance of C v in the IR of our solution has its counterpart in the fact that the superconducting 
state has a non trivial 9^ at the UV boundary. In the next section we investigate this connection in 
more detail, in particular the AdS/CFT picture of our model is completed if we can describe what is the 
role of Ooo in the dual field theory. As we will see, the connotation of 9^ in the dual field theory is the 
final piece of our story. 

4 Dual field theory and marginal deformations 

Having studied the model in the variables {77, 9} it's convenient to go back to the original ones {21, 22}. 
The asymptotic expansion for the two complex scalar fields was given in the (|2.18[) . They both have the 
same near-boundary behavior, 



z l(r ) = ^l_ + ^ + ... , Z2(r) = i^ + ^_ + ... . (4.66) 



We use the definition of z\ and 22 in terms of r\ and 9, i.e. z\ (r) = rcos(0/2) and z%[r) = Tsin(0/2), 
to relate the coefficients #oo, £ and O n to the four parameters appearing in (|4.66[) . The result is the 
following, 

0[ 1] = \0 v caJ-f , 0j» = -\o^sm?f , (4.67) 

0W = i0„am^, 0f= hj^coJ-f. (4.68) 

Given a superconducting solution with fixed 9^ none of the above coefficients is independent of the 
temperature. Nevertheless, defining the constant A = tan(0 oo /2), the following relations hold, 

0« = (4.69) 
Of = -XOi 2) (4.70) 

O^O^ +G ( 2 ) { 2 ) =0 (4.71) 

We observe that A has zero dimension. From the expressions (I4.69[) and (|4.70l) we recognize all the mark- 
ings of the double trace deformation in the AdS/ CFT setup [46] . In the following we explain the details 
of this deformation. They are important in order to understand the dual description which is behind 
the relations (I4.69P and (|4.70[) . A useful technique in this context is the holographic renormalization 
approach. In particular, we want to evaluate the euclidean action on bulk a solution. 

Standard methods allow one to compute the euclidean action Sb as a total derivative (4"5H4"5] . 



S E = - d 



j4, 



~ gC = 2^2 / Sx [°° dr d r (2rh(r)f(r)) = ± I d 3 x (2rM0/(r)) | . ( 4 -72) 



x 



where C is given in (|2.2[) . The surface term at the horizon vanishes both at finite temperature and at 
zero temperature. At finite temperature the term rhh(rh) is finite but f(rh) = by construction. At zero 
temperature the term rh(r) is bounded by C n in the limit r — > but f{r) vanishes like r 2 . The surface 
term at infinity, i.e. in the limit r —> 00, is not finite and needs to be renormalized. At the boundary we 
find two types of divergences, 

(2rh(r)f(r)) = 2r 3 + V r - 2M + | V O^Of . (4.73) 

V / r— 7-00 ^ — ' O — ' 

i=l,2 i=l,2 
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The r 3 term originates from the integration over the AdS space and it is regulated by the Gibbons- 
Hawking term plus a boundary cosmological constant, 



i r* , — 4 



SgH= 2^J d x { 2K+ LJ ' (4 ' 74) 

The metric gs is the induced metric at the boundary and K is the trace of the extrinsic curvature defined 

by 

p^-i(VV + VY) , (4.75) 
with outward pointing unit vector, normal to the boundary. 

The term which is linearly divergent in r comes from the integration over the radial profile of the scalar 
fields. It can only be removed adding boundary counterterms explicitly built out of the scalars z\ and zi- 
They are not unique unless the choice of quantization scheme is fixed. In the present case, the mass of 
the scalar fields is m 2 L 2 = —2 and therefore both scalars can be quantized in the two possible schemes. 
However, there is no ambiguity. Indeed, we recall that £i and C2 are part of an hypermultiplet in the 
TV = 2 theory. Thus, their quantum numbers in AdS are A = 1 and A = 2 and the natural choice is to 
quantize them in different ways. For concreteness we consider the alternative quantization scheme for z\ 
and the standard quantization scheme for zjj. 

The scalar field z\ is dual to an operator of dimension A = 1 whereas the scalar field z^ is dual to 
an operator of dimension A = 2. The boundary values and are interpreted as condensates 
whereas, and are the sources. This choice is in agreement with the dynamics studied in section 
I3.1l for the value 9^ =0. In this case we have set the scalar field £2 to zero and the condensation has been 
associated to the operator with dimension A = 1 dual to £1. According to this choice, the counterterm 
for the Z2 scalar is 

whereas the counterterm for the z\ scalar is 

Si = J d?x^g^(A zm^d^zi + 2 zl/L) . (4.77) 

The renormalized euclidian action is finite and it is given by, 

Sren = S E + S GH + (Si +Si) = ±J d*x + ioPoP - foW) . (4.78) 

Taking into account the relations (I4.69[) and (|4.70[) . the final result is, 



Sren = 4 / <Px ( ^ + 2AO^of ) . (4.79) 



'M 

K 2 J " ~ V 2 

From the above expression we understand the consequence of having identified the sources with the 
condensates. The double trace deformation O^O^, with marginal coupling A, shows up as a finite 
contribution to the renormalized action. In order to be completely general, we note that the identification 
between sources and condensates given in (I4.69P and (|4.70l) takes into account only the "radial" part of 
the complex fields £1 and £2- Indeed, looking for black hole solutions we dropped the phases from our 
Lagrangian however, the double trace deformation has to involve complex operators. We can easily restore 
the phases because the most general solution, according to our ansatz (|3.26p . has the phases ip and <p 
constants. Therefore we simply consider that (j = zi(r)exp(z(<p + ip)/2) and (2 = Z2(r)exp(— i(ip — ip)/2). 
Then, from the asymptotic expansion of Zi(r) and 22(f) given in (|4.66[) . we obtain complex sources and 
complex condensates. The relations (|4.69[) and (|4.70[) are generalized taking into account the following 
observation: it is only consistent to identify operators with the same quantum numbers thus, the charges 
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of the scalar fields fix the relation between sources and condensate. This means that (,1 is related to Q- 
Proceeding with the identification we find 

C>«t = x c o^ ] , = -A c C> 2 2)t , with A c = Ae"^ . (4.80) 

The marginal deformation is then, 

5S = f (PxfXcOxOz + XcO\0^j , (4.81) 

where we have dropped the upper indexes and we have defined {©1,02} as the complex operators dual 
to the scalars {Ci,C2}- Remarkably the deformation is invariant under the U(l) action, 

O x -s- e ia Oi , 2 e- ia 2 , (4.82) 

and this global U(l) symmetry of the field theory is not explicitly broken by 5S. It is worthwhile to 
mention that the operator O^O^ is strictly marginal only in the large N limilH]. 

The understanding of 9^ in the gravitational description is now clear: the parameter 9^ is mapped 
to the marginal coupling A in the dual field theory through the relation A = tan(0 oo /2). It is therefore 
tempting to consider the dependence on the extra coordinate r of the field 9(r) in an RG fashion. In 
fact, it is a general feature in the AdS/CFT scheme to interpret the evolution of the (zero temperature) 
geometry, from the AdS boundary towards the bulk, as the flow of the UV microscopic theory towards a 
low energy regime. In this sense the interpolating solutions represent a novelty of this paper: whenever 
#00 7^ the theory enters a low energy confining phase. However, the double trace deformation that we 
have identified does not break conformal invariance and therefore cannot be responsible for the RG flow. 
On the other hand, it is certainly evident, from the analysis of the potential in Figure [TJ that 9^ 7^ 
drives the theory towards the 9 = n/2 well, away from the conformal fixed point represented by &. This is 
a consequence of 6 being a saddle and cannot be avoided in the classical approximation to supergravity. 
Therefore, we conclude that conformal invariance is somehow broken. To approach this issue we consider 
the asymptotic behavior of 9(r), 

0(r) = 0oo + - + • • ■ . 

r 

In a superconducting solution we have seen that £ = implies 9(r) constant in the bulk. The value 6*00 = 
belongs to this case and even if it represents a trivial example, because A = 0, it matches the expectation 
that conformal invariance is not explicitly broken by the deformation. Therefore, £ ^ parametrizes 
our ignorance of how conformal invariance is broken in the UV |50j . This observation suggests that a 
relevant deformation is turned on at a certain high energy scale. The same relevant deformation is then 
responsible for the RG flow. The intuition on £ ^ can be further specified by considering another 
important relation which holds between the two condensates: 

{ 2 } = . (4.83) 

This is an "on-shell" relation which is not specified by the boundary data given in (I4.69[) and (|4.70[) . 
Instead, it arises from the bulk dynamics of the interpolating solution, in particular from (|4.67p and 
(I4.68p . Because we want to keep the U(l) invariance manifest, it is convenient to work with the complex 
version of (|4.83[) . 

= -CO^e-^ . (4.84) 

According to our interpretation of SS as double trace deformation, the coefficient £ assumes the role of an 
energy scale proportional to the beta function of the coupling A, i.e. r9'{r) oc £/r. We conclude that the 
relation (|4.84l) . which introduces the UV scale £ in the definition of the condensates, breaks conformal 

II We thank Mukund Rangamani for a remark about this point. 
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invariance and provides the relevant deformation which initiates the RG flow. Indeed, by substituting 
the above relation (|4.84p in SS, we find an effective UV action of the form, 

SS ef f oc \^0\O x . (4.85) 

This effective action has the expected RG behavior and in fact, by using dimensional analysis, we find 
that J d 3 x 5S e ff oc (£,/E). Once the RG flow is initiated the coupling A runs. Its value in the IR, given 
by 9(0) = 7r/2, characterizes the effective low energy confining phase with respect to the conformal point 
identified by @. In particular, the field theory dynamics of the tachyonic direction at the saddle point 
can be explored with the following argument. The coordinates (1 and (2 also diagonalize the Hessian in 
a neighborhood of &. The £2 direction corresponds to a scalar with the negative mass m 2 L 2 = — I2/tF*I. 
dual to a relevant operator in the IR conformal theory. Then, any mixing between (1 and C2 in the 
UV will source this IR relevant perturbation. In other words, even if £ is of order e, and the theory 
stays arbitrarily close to the conformal point & at an intermediate scale, the mentioned IR relevant 
perturbation will drive the flow along the tachyonic direction. The above intermediate scale is C v . This 
fact can be seen both from the radial component of the metric in (I3.50P and from the log term that 
appears in the IR expansion of r](r), i.e. rj(r) w \og(C 2 /r 2 ) for r -C C v . 

Finally, the deformation SS is exactly marginal only for the two cases 9^ = and 9^ = it/2. It 
is therefore interesting to note that when 9oo 7^ the coupling flows along the flat direction towards 
9 = tt/2. If 9 = 7r/2 was a second fixed point, this kind of flow would be the expected running of an 
exactly marginal operator when some massive field is integrated out |52) . It is a nice property of our 
model that 9(0) — tt/2 still characterizes the IR effective theory even if it does not correspond to any fixed 
point. In this sense the supergravity interpretation, that 9^, — and 9^ — tt/2 are better understood 
as "order parameters" for the different phases of the model, is also valid in the dual field theory. In the 
next section, we seriously take into account the properties of SS by considering how our model satisfies 
some general statement about Landau- Ginzburg theories. 

Despite the remarkable interpretation of our extremal interpolating solutions in terms of RG flow 
dynamics, a field theoretical argument able to explain the reason why the retrograde condensate does 
not loop back to zero temperature is still absent. It remains an open problem and unfortunately our 
analysis of Section 13.21 only provides other evidence for the non existence of a zero temperature solution 
associated to the retrograde condensate. 

4.1 Thinking about the dual field theory 

We want to suggest a feasible connection between the existence of the moduli space A4g and the existence 
of a marginal operator in the dual field theory. We borrow part of the story from the theory of the N = 2 
Landau-Ginzurg models in two dimensions. Reviewing these ideas we closely follow [53 . 

We recall the notion of moduli in the Landau-Ginzburg model. A pedagogical example, which is also 
useful to understand the nature of phase transitions, is the mean field theory of a single scalar field ip. 
The potential of the theory is given by 

V = m 2 (T)tp 2 + \ip 4 . (4.86) 

with T a continuous parameter (the temperature) and with A constant, positive and greater then zero. In 
a stable configuration the scalar is seated at the minimum of the potential. A smooth phase transition is 
realized if the parameter m 2 (T) takes negative values in some range of temperatures. If m 2 (T) behaves 
as follow, 

m 2 (T) < for T < T c , 

m 2 (T) = for T = T C , (4.87) 

m 2 (T) > for T > T c , 

Note that the mass value is above the BF bound and therefore S is a stable non supersymmetric fixed point 1511 . 
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the extrema of the potential, which is a solution of the equation, 



d v V{T lV ) =0, (4.88) 

changes with the temperature. Above the critical temperature the potential is a sum of positive quantities 
and the only solution to (|4.88l) is <po = 0. Below the critical temperature the value ipo — becomes a 
local maximum and the new minimum, 



^ = V~^' (4 - 89) 

is the stable configuration. We want to emphasize a feature of V(T C ) that is physically crucial. At the 
critical temperature the equation 

d v V{T c , if) = \<pl = (4.90) 

Vo 

has solution ^0=0 which is three times degenerate: we say that V(Tc, f) is critical. Starting from 
this example we analyze a more general situation. The definition of a critical potential falls into the 
framework of singularity functions. By considering n field variables, we define a (polynomial) potential 
V(Xi, . . . , X n ) to be critical, or a singularity function, if its critical points are degenerate. 
The possible deformations of V(ATi, . . . , X n ) are polynomials in the original field variables and are classi- 
fied according to the renormalization group as: relevant, marginal and irrelevant. Relevant perturbations 
split the degeneracy, but marginal deformation do not. The presence of marginal operators reveals that 
the critical potential is not isolated, but rather it is an element of a family My of functions. If k is 
the number of marginal operators then the elements of My are labeled by k continuous parameters 
{Ai, . . . , Afe} which are the coefficients of the marginal perturbations. In other words, an element of My 
is of the form, 

k 

V{X 1 ,...,X n> X 1 ,...,X k )=V(X 1 ,...,X n ) + J2 *i Ti , (4.91) 

i 

where T% is a marginal operator. The space My has the structure of a ring and the coefficients 
{Ai,...,Afc} are called moduli. Each potential living in My admits a critical point which is degen- 
erate. 

If we want to relate the degeneracy of the potentials in My to our theory, Figure |2~T1 clearly shows 
that all the condensates arises from the same branch at T « 0.121, independently of 6*00. We can say 
more about the relation between the spaces My and Me- Taking into account the results of the previous 
section we can collect the following chain of observations. First, the double trace operator T = O1O2 
is a marginal operator which belongs to My in the dual field theory. Second, the modulus A associated 
to T is geometrically the parameter which belongs to Me- Thus, the AdS/CFT correspondence 
maps the modulus associated to JF to the space Me- More precisely, the relation Ac = tan(0 co /2)e -1 ^, 
together with the restriction Ooo € [0,7r/2], which defines Me, implies that Ac parametrizes the unit ball 
on the complex plane. We remove from this set the circle S 1 because of the retrograde condensate and 
we refer to the open ball as B(0, 1). At this point we can rephrase the analysis of the previous section by 
considering the action of the renormalization group flow on the open ball. We know that the origin is a 
fixed point and therefore we focus on the set B(0, 1) \ {0}. In this case the renormalization group flow 
acts as an projection and the IR image of £>(0, 1) \ {0} is the unit circle. 

At the classical level one might expect Ac to vary in the entire complex plane. Instead, we find that 
Ac G B(0, 1). It is natural to ask what happens to the complement of B(Q, 1). In this case 800 takes 
values in the range [tt/2, it]. We already know the dynamics of the model simply because the potential is 
7r/2-periodic: when 6^ is increased from 7r/2 to 7r, the condensate goes backwards from the retrograde 
condensate to the conformal domain wall. These solutions are related to the ones found for 6^ e Me but 
they are not the same. In fact, for 9^ = it the scalar field £1 is set to zero and the condensation is driven 



23 



by (,2- This situation is the opposite of the case 9oo = 0. Let's see what is the boundary description in 
this case. The relations 

= XO { l ] , Of ] = -XO^ , (4.92) 

are always valid, even in the range € [tt/2,tt]. In the limit 9oo — > tt, the coupling A blows up and 
the relations (|4.92l) make sense only if (D± — > and — > 0. These two conditions are suitable for 
the opposite quantization scheme from that adopted in ,8(0, 1): the condensates are now and 0f\ 
the sources instead are and . Thus, the conformal domain wall in the case 9^ — -k has to 
be associated with the condensation of the operator O^' with ojp — 0. In summary, the quantization 
schemes in the region C \ B(0, 1) is opposite from that in B(0, 1). 



5 Conclusions and Outlook 

In this paper we have studied a particular N = 2 supergravity theory in four dimensions. We based our 
analysis on the general principles of the AdS/CFT correspondence focusing our attention on the holo- 
graphic superconductivity framework. The theory has interesting dynamics and several new ingredients 
coexist in the final picture. These are, the interpolating solutions constructed in section 2, the properties 
of the zero temperature solutions outlined in section 3 and the notion of double trace deformation in 
AdS/CFT. Each of them can be related to the existence of the space Aig which is completely tied to the 
nature of the hypermatter scalar manifold SU(2, l)/t/(2). We emphasize the topological properties of the 
manifold SU(2, 1)/U(2) in the construction of the model. The manifold SU(2, 1)/U(2) is homeomorphic 
to a ball in C 2 and the coordinates (|3.26p . for fixed r / 0, are indeed related to the Hopf fibration of 
the three-sphere. In this description the angle 9 is the azimuthal angle of the two-sphere whereas the 
polar angle if and the fiber ip are phases. Phases cannot appear in the potential and the potential V is 
function only of r and 9. At the origin, r = 0, the compact space shrinks. This fact manifests itself in a 
"topological" degeneracy and Me shows up as a moduli space. 

M = 2 supergravity puts strong constraints on the nature of the scalar manifold: hypermultiplets 
parametrize quaternionic Kahler manifolds and vector multiplets parametrize special Kahler manifolds. 
For the quaternionic case, the most relevant manifolds are the homogeneous spaces of real dimension 
4m [51ES], 

„. . SU(2,m) , „„. 

X{m) = SU(m) x SU{2)x U(l) " ^ 

The lowest dimensional space corresponds to m — 1 and it is the coset manifold SU(2, 1)/U(2). In 
this sense we can say that the hypermultiplet {(1^2}, considered as building block of our model, is 
universal. Then, the graviphoton together with the metric connections, are the minimal elements needed 
to construct the Lagrangian (|2.2[) . It may be possible that the topology of higher dimensional quaternionic 
manifold leads to a more general moduli space. In this case the gauging of the U(l) isometries requires 
the presence of vector multiplets. In four dimensions vector multiplets contain complex scalars and the 
interplay among all the matter fields can be quite involved. If the deformations T% £ M.y are constrained 
by some U(l) global invariance, then the relevant gauging can only pick the 03 direction in the SU(2) 
part of the isotropy group. For that reason, candidates to describe a more general moduli space are 
scalar manifolds which contain copies of SU (2) in their isotropy group. On the other hand, the presence 
of several gauge fields in the gravitational action introduces more than one U(l) charge in the dual field 
theory. The complexity of the bulk system is therefore justified by the complicated dynamics of the field 
theory. 

The class of 5D JV = 2 supergravity theories coupled to matter fields is also very well understood. In 
this context, the gauging of the universal hypermultiplet follows the same procedure outlined in section 
[2j When the U(l) killing vector is taken to be the 173 direction in the SU{2) part of the isotropy group, 
the two complex scalar fields parametrizing the quaternionic manifold have opposite charges and the 
same mass. In particular, the mass lies at the conformal value m 2 L 2 — —15/4. The model has been 
proposed as a truncation of the type IIB theory on AdS$ x T 1,1 and the two scalars might describe the 
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dynamics of the chiral operators Oij = Tr[AiBj] in the dual field theory [5B]. In the next paragraph 
we make a more precise comment about the implications of having a brane construction of a generic 
model. Concerning this 5D dimensional truncation, we only mention that the model contains a standard 
holographic superconductor and a retrograde condensate [57] ■ In principle, the analysis carried out in the 
present paper has a straightforward generalization. Indeed, an extremal solution similar to (|3 . 50[) exists 
also in this five dimensional case and is given by, 



However, the radial component of the metric, g rr = (/(r)/i 2 (r)) _1 , goes to zero like r in the limit r — > 
and therefore the IR geometry is not a cone. It is not clear whether this behavior could affect substantially 
the physics of the model. 

The Lagrangian (|2.2p belongs to the category of the "top-down" models. With this name it has 
been indicated the set of all those theories that can be constructed by first principles. In general, a 
top-down model is a model which admits a description in terms of branes, either in the type IIB theory 
or in the M-theory. Having a brane description of the theory allows a more precise understanding of the 
relation between bulk fields and boundary operator. This statement has a strong physical implication, 
it means that the theory has a microscopic description in terms of open strings. Thus, the potential 
and other couplings of the model are fixed by the string dynamics. Several examples have been found 
and discussed [58ft6T] . In our case the Lagrangian (|2.2I) can be embedded in a consistent truncation 
of the four dimensional Af = 8 supergravity however, the relation with the ABJM theory encounters a 
spontaneous obstacle. Indeed, a remark in £37: shows that the boundary conditions for a U(l) spontaneous 
symmetry breaking superconductor, allowed by this theory, are restricted. In particular, if we insist on the 
spontaneous symmetry breaking, the microscopic theory for our boundary conditions (|4.69p and (|4.70p . is 
not known. The present paper provides evidence of a beautiful connection between the scalar manifold of 
the gravitational theory and the space of marginal deformations of the dual field theory. If the dual field 
theory admits a brane construction it might be possible to provide a general statement about this type 
of connection. In this case the supergravity Lagrangian is obtained by KK reduction and the scalars of 
model parametrize the complex structure of the compact internal space, on the other hand this compact 
internal space has to do with the moduli space of the field theory. These kind of ideas deserve further 
investigation. 

Finally, we want to stress the importance of the retrograde condensate. It lives at the border of 
the moduli space and even if it is not relevant to the thermodynamics of the system, it does signal 
the existence of more general solutions. Furthermore 8 = tt/2 is also the IR limit of the renormalization 
group flow of the marginal coupling constant and the value for which the charged soliton and the extremal 
solution (|3.50p have been found. Another example of retrograde condensate arises in the context of Af = 8 
supergravity in five dimensions whose dual field theory is Af = 4 SYM [62]. The sector of the Af = 8 
theory in which the retrograde condensate has been found is the Af = 2 Lagrangian of the STU model 
coupled to three hypermultiplets [53] . We suggest that the theory could have an interesting structure not 
yet explored. The study of this case would be particularly relevant for understanding the phase diagram 
of Af = 4 SYM. 
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$(r) = , 6{r) = - . (5.94) 
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A Black hole solutions: horizon expansion and thermodynamics 



The superconducting black holes are obtained integrating the equation of motions from a postulated 
horizon. It is a numerical problem and in order to be well posed we have to discuss how the field 
variables are expanded around this horizon. We use the variables {f(r),h(r),$(r),rj(r),9(r)} that are 
the ones in which the black hole solutions have been constructed (see section l3~Tj) . The general strategy 
to attack the problem is explained in |13j . here we repeat the main ingredients. 

A regular event horizon is defined by the greatest simple root of the metric components, we call this 
root r/j. The temperature of the black hole is the Hawking temperature, 



T = 



i- yjf'(r)(f(r)h*(r)y 



(A.95) 



The metric component gu and g rr have as a common zero. For the numerical problem we can choose 
the gauge 

f(r) = f(r h )(r ~ r h ) + . . . , h(r) = h(r h ) + h'(r h )(r -r h ) + ... . (A.96) 

in which f(rh) = but h(rh) is finite. Being the horizon a simple root we need to consider /'(r^) 7^ 0. 
The equation of motion for the Maxwell field $ involves the time component of the contravariant quantity 
A 1 - 1 . For A 1 to be well defined we require $(r) to vanish like f(r), 



$(r) = $'(r h )(r - r h ) + 



(A.97) 



Then the equation (|3.33[) is well defined at the horizon. The two scalar fields, 77(f) and 9{r), are assumed 
to be constant at the horizon, therefore their expansion is 



r](r) = rj(r h ) + rf '{r h )(r - r h ) + 



{r)=9(r h )+9'(r h )(r-r h ) + .. 



(A.98) 



In this case the equation of motion (|3.34|) and (|3.35[) are well defined only if the divergent term, propor- 
tional to 1/ f[r), cancels. This requirement imposes a condition on the value of the first derivatives of 
the scalars field. In this case, rf{rh) and 9'(rh) are given by, 



1 



d V 



sinh 2 (77/2) f'h 2 



(A.99) 



The coefficient f'(rh) is not an independent quantity, in fact its value is fixed by the consistency of the 
equation (|3.32p . 



2 



(A.100) 



Finding the solution to the equations of motion is now a well posed Cauchy problem. We can list how 
many free parameter are left, 



h(r h ) , $'(r fc ) , 77M , 9(r h ) . 



(A.101) 



This set is further reduced by considering that the equations of motion are invariant under the two scaling 
symmetries, 



t -> at , <f> -> $/a 



/ f/a 2 , 



ah 



(A.102) 



r^ar, {t,x,y) -> (t,x,y)/a , f-+a 2 f, $ -> a$ . (A.103) 

The scaling (|A.102p is used to set h(rh) = 1 and the scaling (|A.103[) is used to fix the value of the horizon 
coordinate, we take rh — 1. Then, we have three independent parameters, {^(r^), T)(rh), 9(rh)}, and only 
two asymptotic conditions to match. It follows that the solution is specified only by one parameter. Once 
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the charge of the black hole has been fixed to a constant value, for example p = 1, this parameter is the 
temperature (|A.95|) . 

The superconducting black holes are dual to a field theory ground state. The thermodynamical 
properties of this ground state can be investigated using the bulk solutions. We briefly review how the 
field theory thermodynamics is recovered from the AdS/CFT correspondence. 

The entropy density s of the boundary field theory is related the surface area of the horizon through 
the formula, 

a = ^ . (A. 104) 

The other thermodynamical quantities of interest are: the energy density e, the free energy /, the charge 
density p and the chemical potential p. They are related to the coefficients of the boundary expansion, 
(|2.15p and (|2.f 7p . and are given by, 



M „ . _ pp 

2k 



PP=^, f = e-Ts. (A.105) 



In the canonical ensemble the following thermodynamical identity holds, 

/ = e- Ts = Cl + pp , (A.106) 

where f) is the (density of) grand canonical potential. This potential is calculated from the action S ren 
which has been introduced in Section^] We can directly substitute into (|4.73p the values of explicitly 
given in (|4.67[) and (|4.68p . where i, j — 1,2. This step yields the result, 

(2rh(r)f(rj) = 2r 3 + -Ol r - 2M . (A.107) 

V / r— >oo 4 

Therefore, in order to rcnormalize the action we need to add the Gibbons-Hawking term plus a boundary 
cosmological constant and the countcrtcrm for a A = 1 condensate dual to rj(r), 

Si = --^5 J dW=ffi( ^K^) + r, 2 /2) . (A.108) 

The final result gives the relation -S ren = J d 3 x M/(2n 2 ) and thus, Cl = TS E = -M/2. 
The relation (|A.106I) can be derived from a conserved charge. The procedure to obtain the conserved 
charge is straightforward once the Einstein equation, R xx — K 2 g xx C = 0, is considered. This equation is 
not independent from the set of equations (|3.31l) - (|3.35l) and can be manipulated to prove that, 

Q(r) = r 2 h(r)f(r) - 2rh{r)f{r) - r 2 /i(r)$(r)$'(r) (A.109) 

is a constant of motion of each black hole solution [34]. The quantity Q corresponds to the Noether 
charge of the scaling symmetry (IA. 103|) 36! . Then, the relation Q(r^) = Q(r —> oo) is equivalent to the 
thermodynamical identity (|A. 1061) . 



B The Superpotential 

We derive the solution (|3.50p from a superpotential equation. This is general procedure that has a 
natural derivation in supersymmetric theories. In the following we review the steps of this construction. 
We proceed as in [64) . setting $(r) = and f(r) = r 2 . Under certain assumptions, the main idea is simply 
to arrange the second order equations of motion in a set of first order equations. One of the equations 
introduces a new quantity W called superpotential. The important point is that the superpotential is 
defined by an equation which involves only the potential V and its derivatives, d^V and dgV. 
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We are interested in the equations (|3.32p . (|3.34p and (|3.35p . Defining the field variable g(r) = 
f(r)h(r) 2 — r 2 h{r) 2 the equations take the following form, 

V-isinh 2 (2)^ + 4+1 
4 4 \2J r 2 2g 



- A V 12 - - smh 2 ( i ) 0' 2 + - + — Pfa, 0) - , (B.110) 



n" + ( - + 7t) ff ~ 7 si nh »? #' 2 - -fyPfa, 9) = , (B.lll) 
\r 2g) 4 cj 



(B.112) 



+ — ) 9' + coth(^)n'6' i -d e V(-n,0) = . 

2gJ \2J 1 sinh 2 (77/2)3 * Kh > 

The first one is a constraint and can be solved for g(r), 

g(r) = P( V , 9) Qr/ 2 + ~ sinh 2 (|) tf' 2 - -|) . (B.f 13) 



Then, substituting the above expression in the remaining equations, (jB.llip and (jB.llip . we obtain two 
equations for the two field variables, r\(r) and 9(r). The final result is, 



n 



\ - T - sinh 2 (|) 9' 2 \ rf - - \ sinhr, 9' 2 + ?£z = 0, (B.114) 



1 - V + coth 2 ("I) A 9' T - sinh 2 (2) 9'" + * = , (B.115) 

r 4 ' V2/ 'J 4 V2/ sinh 2 (r?/2) V y ' 



where, 



At this point, the second order equations of motion can be reduced to first order equations. This procedure 

is general and can be written in a formal way: given the superpotential W the set of three non linear 
first order equations iJ^l. 

,'=-^, (B.H7) 
r Vv 

9>=- A - I (B.118) 

r sinh 2 (77/2) W 7 

7> = 4(<9 I; W) 2 + 2 4 - (d e W) 2 - 3W 2 . (B.119) 

sinh (?7/2) 

In our model the superpotential can be found by inspection of the equation (|B. 1 19|) or directly from the 
supersymmetry transformations of the fermions [26] . It is given by, 



W = \j 2 cosh 2 (|) + ~ cos 2 9 sinh 4 (|) . (B.120) 
Then, it is immediate to check that 



C 7T 

n(r) = 2 arcsinh , 9(r) = - , (B.121) 
r 2 



is a solution of the system (|B.117HB.118p . 

ttTo match the notation of [64] we look at their Lagrangian (2.1). Then r t -> 8 -> ^28 with V ^ IV and W -> v^W 
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